Optimal time-dependent polarized current pattern for fast domain wall propagation in 
nanowires: Exact solutions for biaxial and uniaxial anisotropies 
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One of the important issues in nanomagnetism is to lower the current needed for a technologically 
useful domain wall (DW) propagation speed. Based on the modified Landau-Lifshitz-Gilbert (LLG) 
equation with both Slonczewski spin-transfer torque and the field-like torque, we derive the optimal 
spin current pattern for fast DW propagation along nanowires. Under such conditions, the DW 
velocity in biaxial wires can be enhanced as much as ten times compared to the velocities achieved 
in experiments so far. Moreover, the fast variation of spin polarization can help DW depinning. 
Possible experimental realizations are discussed. 
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Introduction.— Fast magnetic domain wall (DW) prop- 
agation along nanowires by means of electrical currents 
is presently under intensive study in nano-magnetism, 
both experimentally [H-Q and theoretically In ad- 

dition to the technological interest such as race track 
memory DW dynamics is also an interesting fun- 
damental problem. The dynamics of a single DW can 
be qualitatively understood from one-dimensional (ID) 
analytical models |9l-lll| that predict a rigid-body prop- 
agation below the Walker breakdown and an oscillatory 
motion above it 0, 12]. The latter process is connected 
with a series of complicated cyclic transformations of DW 
structure and a drastic reduction of the average DW ve- 
locity. The Walker limit is thus the maximum velocity at 
which DW can propagate in magnetic nanowires without 
changing its inner structure. Prom a technological point 
of view, such a limit seems to represent a major obsta- 
cle since the fidelity of data transmission may depend on 
preserving the DW structure while the utility requires 
speeding up the DW velocity adequately. Various efforts 
have been made to overcome this limit through geome- 
try design. For instance, Lewis et al. [1 31 ] proposed a 
"chirality filter" consisting of a cross-shaped trap to pre- 
serve the DW structure. Yan et al. [14J demonstrated 
the removal of Walker limit via a micromagnetic study 
on the current-induced DW motion in cylindrical Permal- 
loy nanowires. In this Letter we investigate other ways 
to substantially increase the DW velocity avoiding the 
Walker breakdown. 

A DW propagates under a spin-polarized current 
through angular momentum transfer from conduction 
electrons to the local mag netization, known as the spin 
transfer torque (STT) [15| . which is different from mag- 
netic field driven DW propagation originated from the 
energy dissipation [l2|, Us]. Generally, two types of 
spin torques are considered: the Slonczewski torque 15 1 
(a-term) T a = -7-^-M x (M x s) and the field-like 
H [3 (6-term) T b = -76./M x s, where 
/m e , M, M s — |M|, and s are the the gyro- 



torque 
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ration magnetization, and the spin polarization direction 
of itinerant electrons, respectively, aj and bj depend on 
current density j e and spin polarization P. Theory pre- 
dicts [HEll that aj = Pj e h/2d\e\M s and bj = f3aj, 
where d is the thickness of the free magnetic layer. /3 is 
a small dimensionless parameter that describes the rel- 
ative strength of the field-like torque to the Slonczewski 
torque. The value of /3 is sensitive to the thickness of the 
free layer and the decay length of the transverse com- 
ponent of the spin accumulation inside the free layer as 
discussed in Ref. I l8l| . The typical value of /3 ranges 
from to 0.5 18j, [l9(. In the conventional case of current 
along the nanowire with biaxial magnetic anisotropy, the 
a— term is incapable of generating a sustained DW mo- 
tion, except for a very large current, while the b— term 
can drive a DW to propagate Q. Unfortunately, the 
b— term is usually much smaller than a— term 0, A 
large current density is needed in order to reach a techno- 
logically useful DW propagation velocity [l| , but the as- 
sociated Joule heating and DW structure collapse could 
affect device performance. We show that the problem 
can be solved if one uses an optimal polarized current 
pattern. 

In this Letter, our focus is on the optimal spin- 
polarized electric current pattern for fast DW propaga- 
tion along nanowires. For usual magnetic materials, our 
theoretical results show that the DW velocity can be en- 
hanced by as large as ten times in comparison with DW 
velocity driven by the conventional constant current in 
existing experiments. Moreover, the ultrafast change of 
spin polarization can be used to de-pin a DW. 

Model.— The internal magnetic energy of a nanowire 
can be formulated as 



U[M] = 



d 3 s 
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sin 



magnetic ratio, magnetization of the magnet, the satu- 
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where 9 and 4> are the polar angle and azimuthal angle 
of the local magnetization m = M/M s , respectively. J 
and w are the exchange energy constant and energy den- 
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sity due to all kinds of anisotropics, respectively. The 
dynamics of M is governed by the LLG equation [2(| : 



— = - 7 M x H e// 



— M x — — 

M s dt 



-STT, 



(2) 



here H e ff = ~j^$U/8bA is the effective magnetic field 
and a is the phenomenological Gilbert damping constant 
[20j| . T S tt is the STT with both the Slonczewski-type 
and the field-like terms. 

Biaxial anisotropy .— Considering a biaxial anisotropy 
w(6,4>) = ~^jm 2 z + 4pm 2 with the easy axis 
along z direction and the hard axis along x direc- 
tion, the effective field takes a form of H e ff = 
j^(Km z z - K ± m x x) + Here K and K ± 

describe energetic anisotropies along easy z axis and hard 
x axis, respectively. We assume that all local spins lie in 
a fixed plane called DW plane, i.e., (z, t) — <f> (t) , which 
should be checked self-consistently late. In the spherical 
coordinates, Eq. ^ becomes 
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where s r , sg, and s^ are three components of unit spin 
vector s in spherical coordinates. The DW profile sat- 



isfies Jfjf - 2»l2£ 
condition of 9 



[K + K± cos 2 0) = with boundary 



and 7r at distance. One obtains the fa- 
mous Walker's DW motion profile tan | = exp (^ z ~ x ^ 



2 ~~ \ A 

in which X (t) is the position of the DW center and 
A = \J k+k± cos 2 <t> ^ s DW width resulting from the bal- 
ance of anisotropy energy and exchange energy Q . These 
assumptions are valid under sufficiently low current den- 
sity which will be demonstrated later. Substituting this 
DW profile into Eqs. ^ and (g]), we have 

A smO 2/iqMs 



x 1 
a— + = 7 



Q j ( s - fog) 
sin 9 
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For a given DW motion, term — /3sg should be as 
large as possible in order to lower the needed current 
density. Meanwhile, considering identity (s^ — (3sg) 2 + 
(s + fis^) 2 = (l + P 2 ) (l - s 2 ) , we choose 



s r = 0, sg = cos 77, 



sin 77, 



(7) 



with optimization parameter rj. 

To ensure the spatial-independence of X and 0, the 
above equations require aj to be proportional to sin9, so 



we let a j = Aj sin 9 — Aj scch( z ) with a constant 
A 7. Thus we have 



AT = 7 A 

= 7 
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which describe the DW propagation and the DW plane 



precession. 



Here 



(V) 



b'j ( 7 ?) = Aj (cos-q + /3 sin 77). 



Aj (sin r) — /3 cos ?/) and 
The DW width is time- 
independent when the DW undergoes a rigid-body propa- 
gation with (t) — 0o =constant. In general, DW width 
A depends on the time through time-dependence of 0. 
The exact rigid-body solutions constitute 
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2^o M, 



sin 20o (77) 
A (57) 




-ff + ld. cos 2 0o (??) 



7— A (r?) . 



(12) 



The spin current pattern is then described by 77. Differ- 
ent value leads to different canted angle, DW width and 
propagation velocity. It is straightforward to show that 
the assumption of rigid-body motion always holds under 
condition Ajy/{l+a 2 ) (l+p 2 ) < 

Before finding the optimized spin current pattern for 
maximal velocity, let's first consider two special cases. 
The conventional case in the existing experiments [21. 22| 
is the constant current density with electron spin polar- 
ization along z— axis, i.e., sg = — sin# and = 0, which 
gives the velocity u\ — l^f- A (n) . It again shows that 
Slonczewski torque is incapable of generating sustained 
DW propagation while the field-like torque can. But the 
velocity is rather small since /3 <C 1 in usual materials. 
However, the DW velocity can be greatly enhanced if 
a— term is involved. This is the case of 77 = It gives 
the velocity U2 — 7 4rA (§) • in typical materials (l9j . 
/3 ~ 0.1, so the velocity is 10 times lager than u\. One 
can see that DW propagation velocity is greatly enhanced 
under a modification of the spin polarization and locally 
minimized current density pattern. 

The maximal velocity AT max = X (77*) at the optimal 
parameter 77* can be found through exact numerical cal- 
culations although a closed analytic form is difficult to 
obtain due to the complexity of Eqs. (fT0|) . (fTTj) . and (fT2l) . 
Factor A = X max / Ul = sin n ' ~g cos ^ measures the 
velocity enhancement. Fig. la is the /3 _1 dependence of 
A for various damping coefficients and typical magnetic 
parameters. It is approximately linear, and insensitive 
to damping parameter a. Fig. lb is the plot of spatial 
distribution of s x , s yi s z and aj for the optimized spin 



current pattern around DW center. We note that s x , s 



3 



and s z vary only near the DW center, and reach fixed 
values away from the DW. A large perpendicular compo- 
nent s y is required to achieve large DW velocity The rea- 
son is that perpendicular spin component induces a large 
effective field H a — -f/-M x s. Thus the DW moves un- 
der the Slonczewski torque with a large component along 
wire axis. This finding is consistent with recent micro- 
magnetic simulations 23[ showing that the DW velocity 
can be greatly increased by applying perpendicular spin 
polarizations. It is also very interesting that locally min- 
imized current density aj is finite only near DW center 
while it becomes zero at distance, which should greatly 
lower the energy consumption. 
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FIG. 1: (Color online) (a) DW velocity enhancement factor 
A versus /3 -1 at different damping coefficients a = 0.01,0.02 
and 0.5. (b) The spatial distribution of x, y, z components 
of the optimal spin polarization pattern and current density 
pattern for maximal DW velocity under a — 0.01 and /3 = 0.1. 
The other parameters are using the materials parameters of 
Permalloy: M s = 8.6 x 10 s A/m, J = 1.3 x 10" 11 J/m, 
K = 500 J/m 3 , K ± = 8 x 10 5 J/mffJ, and a reasonable 
value A j = 25 Oe according to Ref. [23| . 



Depinning.— Besides the advantage of markedly speed- 
ing up the DW velocity, our time-dependent spin current 
pattern also implies a possible way to improve the ef- 
ficiency of DW motion against the pinning effect. The 
argument is attributed to an additional force on the wall 
due to a fast changing spin direction. It is convenient 
for us to treat DW as a quasiparticle [HI with mass 

2SnlM~ 



in 
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pinning. 

pinning force -F p i n 



25l [26| when we deal with the effect of 
Here S is the cross section of the wire. The 
is expressed by the pinning 



potential E and the position X of the wall. Thus for 
small <j>, Eqs. flSJ and © can be simply decoupled and 
result in 



F 



= x — aiK± x 1 dE 

(1 + a 2 ) HqM s m w dX 
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where the temporal variation of DW width is neglected. 
The contributions of current-induced acceleration are 



the last two terms. In usual setup (s along the wire 
axis), the depinning acceleration due to STT is X sp — 

— A 1 ^J- [26|. Thus, one can observe that our opti- 
mal spin current pattern provides a higher depinning ac- 
celeration by 10 times for the current density part since 
parameter f3 is typically around 0.1 19]. The force F 
on the wall does not only depend on the current density 
but also on the time derivative of spin direction. The 
switching-time dependence of current-induced depinning 
follows the last term in Eq. (jT3"]) . For a fast changing 
current, i.e., spin variation rate r <C ^°kI ~ 10~ 12 s, 
the contribution of the time derivative term will be sig- 
nificant. The novelty and importance of our proposal 
are embodied in the nature of ultrafast switching-time of 
spin degree of freedom. It makes a main difference from 
the DW depinning due to current density's rise time (26j 
which is limited by the intrinsic response time of circuit 
27|. 

Uniaxial anisotropy .—Let K± = 0, Eqs. ([5]) and © 
result in 



X 2 



~f 2 A 2 A 2 
1 + a 2 



1 + /3 2 
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with DW width parameter A = yj J/K. 

Thus, the largest possible DW propagation veloc- 
ity is X max = 1 A J jn + P 2 )/ (l + a 2 )A. The conven- 
tional polarization 21, [22[ gives the DW velocity u = 
j A j A. As a result, the DW velocity is enhanced 



under the optimal spin current pattern by a factor of 

/« = \A+(i*& " 



x 



We note that the enhance- 



ment is not so large in the uniaxial wire since both a and 
{3 are far less than 1 in usual magnetic materials [28[ . The 
physical reason lies in that a— term is capable of gener- 
ating a sustained DW motion in uniaxial wire, which is 
different from biaxial case. 

Discussion.— Although the optimal spin current pat- 
tern for maximum DW velocity is found, it is still an 
experimental challenge to generate a temporally and 
spatially varying spin polarized current. Interestingly 
enough, a very recent experiment used spin-polarized 
current perpendicular to a nanowire to manipulate DW 



motion 



22 



There are now at least two types of current 
patterns realizable. The hope is that our capable exper- 
imentalists can one day generate any designed current 
pattern. Indeed, there are many theoretical proposals 
for generating a designed current pattern. Tao et al. [29[ 
and Delgado et al. 'AO have proposed to use magnetic 
scanning tunneling microscopic (STM) tip above a mag- 
netic nanowire to produce localized spin-polarized cur- 
rent. Experimentally, the control of spin-polarized cur- 
rent in a STM by single-atom transfer was demonstrated 
very recently by Ziegler et al 3l|. In summary, our pro- 
posed optimal spin current patterns are difficult to gen- 
erate now, but their existence does not violate any funda- 
mental laws and principles. Our results and calculations 
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33j, the [12] 



will be relevant to experiments when the generation of an 
arbitrary spin-polarized current pattern becomes true. 

In the above discussions, the spin pumping effect on 
the DW motion is neglected because the DW-motion in- 
duced current is zero in biaxial wire since there is no DW 
plane precession (32j below Walker breakdown and it is 
much smaller than the applied external spin-polarized 
current in uniaxial wire. According to Ref. 
maximum DW-motion generated electric current density 
in uniaxial wire is (j 2 ) = \ (<j-|- — crjj x ^ ; where L is 
the length of the nanowire, and o± denote the conduc- 
tivities of the majority and minority electrons. In the ex- 
periments of Beach et al. L ~ 40 /im, A ~ 20 nm, and 
DW velocity X max ~ 40 — 100 m/s. For a typical conduc- 
tivity cr t ~ 10 6 tt- 1 m , one can find the pumped elec- 
tric current density no more than (j z ) ~ 10 5 — 10 6 A/m 2 , 
which is much smaller than the typically applied current 
of the order of 10 10 — 10 12 A/m 2 in experiments 

Conclusion.— We propose an optimal spin current pat- 
tern for high DW propagation velocity in magnetic 
nanowires. In uniaxial wires this enhancement is of mod- 
est size, while in biaxial wires a factor of a few tens can 
be achieved. The nature of ultrafast switching-time of 
spin degree of freedom proves to be a novel way to im- 
prove the efficiency of DW motion against the pinning. 
We expect our proposal will stimulate and also possibly 
guide future experiments. 
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